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Abstract. We study the eigenvalue problem — u" + V(z)u = Xu in the complex plane with 
the boundary condition that u(z) decays to zero as z tends to infinity along the two rays 
argz = — ~ ±^L_, where V(z) = —(iz) m —P(iz) for complex-valued polynomials P of degree 
at most m — 1 > 2. We provide an asymptotic formula for eigenvalues and a necessary and 
sufficient condition for the anharmonic oscillator to have infinitely many real eigenvalues. 
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1 Introduction 

Theory of the so-called PT-symmetric Hamiltonians has been a very active research area in 
recent years [H EJ El EJ El El QUI E3]. While they are not self-adjoint in general, many PT- 
symmetric Hamiltonians have real spectra [H El El E]- Self-adjointness of Hamiltonians is not 
physical requirement, but guarantees reality of spectrum and conservation of probability. So 
one searches for non-self-adjoint Hamiltonians that have real spectra. In this paper, we study 
the Schrodinger eigenvalue problems in the complex plane with complex polynomial potentials 
under decaying boundary conditions along two fixed rays to infinity. We will prove that the 
Schrodinger eigenvalue problem has infinitely many real eigenvalues if and only if the eigenvalue 
problem is 'PT-symmetric or it is a translation of a PT-symmetric problem. 
For an integer m > 3 we consider the Schrodinger eigenvalue problem 



Hu(z) :-- 



'dz 2 



P{iz) 



u(z) = Xu(z), for some A 6 



(1) 



7T 2-7T 

2 m + 2' 



(2) 



with the boundary condition that 

u(z) — >• exponentially, as z — > oo along the two rays arg(z) = 

where P is a polynomial of degree at most m — 1 of the form 

P(z) = aiz m ~ l + a 2 z m ~ 2 + ■ ■ ■ + am, a:= (oi, a 2 , . . . , a m ) G C m . 

The existence of infinitely many eigenvalues is known. Sibuya [12] showed that the eigenvalues 
of H are zeros of an entire function of order p := | + ^ € (0, 1) and hence, by the Hadamard 
factorization theorem (see, e.g., [1]), there are infinitely many eigenvalues. 



*This paper is a contribution to the Proceedings of the 5-th Microconference "Analytic and Algebraic Me- 
thods V". The full collection is available at |http://www.emis.de/journals/SIGMA/Pra guc2009. html| 
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The Hamiltonian H in ([JJ with the potential V(z) under the boundary condition ([2]) is called 
PT -symmetric if V(—z) = V(z), z £ C. Note that V{z) = —{iz) m — P(iz) is a PT-symmetric 
potential if and only if a £ W 71 . 

In this paper, we will derive the following asymptotic expansion of the eigenvalues A n , ordered 
in the order of their magnitude. 

Theorem 1. There exists N = N{m) £ Z such that the eigenvalues {A n }^? =Ar of H satisfy that 



where cq{cl) = it 1 sin (^) B (|, 1 + ^) and c\(a) = and Cj(a), 2 < j < m + 1, are explicit 
non-zero polynomials in the coefficients a of the polynomial potential, defined in ()18p . 

As a consequence of ([3]), we obtain a necessary and sufficient condition for the anharmonic 
oscillator H to have infinitely many real eigenvalues as follows. 

Theorem 2. The anharmonic oscillator H with a potential V{z) has infinitely many real eigen- 
values if and only if H with the potential V(z — zq) is VT -symmetric for some zq £ C. 
Moreover, no H has infinitely many real and infinitely many non-real eigenvalues. 

If H is PT-symmetric, then u(z) is an eigenfunction associated with the eigenvalue A if 
and only if u(— z) is an eigenfunction associated with the eigenvalue A. Thus, the eigenvalues of 
a 'PT-symmetric Hamiltonian H either are real or come in complex conjugates. Also, magnitude 
of the large eigenvalues is strictly increasing [10]. Thus, if H is PT-symmetric, then eigenvalues 
are all real with at most finitely many exceptions. 

This paper is organized as follows. In Section [21 we will introduce results of Hille [6], 
Sibuya [12], and Shin [101 [TT] on properties of the solutions of ([I]). We will present asymp- 
totics of solutions of (pQ), introduce the spectral determinant of H, and provide an asymptotic 
expansion for the spectral determinant in the sector where all but finitely many eigenvalues lie. 
In Section [3l we will prove Theorem Q] where we use the asymptotic expansion of the spectral 
determinant. In Section HI we will prove Theorem [2] and provide other spectral results; expres- 
sing A n as an asymptotic formula in n and obtaining an asymptotic formula for nearest neighbor 
spacing of eigenvalues, monotonicity of their magnitude, and more. 

2 Properties of the solutions 

In this section, we will introduce results of Hille [6], Sibuya [12], and Shin [10\ [TT] on the 
properties of solutions of ([I]). 

First, let u be a solution of (pQ) and let v(z) = u{—iz). Then v solves 



Solutions of have simple asymptotic behavior near infinity in the complex plane and in order 
to describe the asymptotic behavior, we will use the Stokes sectors. The Stokes sectors Sk of 
the equation ^ are 




(3) 



j=0 



v"(z) + [z m + P{z) + X]v{z) = 0. 



(4) 




Hille O § 7.4] showed that every nonconstant solution of either decays to zero or blows 
up exponentially in each Stokes sector. Moreover, if a nonconstant solution v of flU decays in 
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a Stokes sector, it must blow up in its two adjacent Stokes sectors. However, when v blows up 
in a Stokes sector, v need not be decaying in its adjacent Stokes sectors. Thus, the decaying 
boundary condition ([2]) becomes that v(z) — > as z — > oo in the sectors Si U S_i. 

We will also use the following functions of the coefficients a € C m of the polynomial P. Let 

b jlk (a) be the coefficient of z mk ~ j in (A (P(z)) k for 1 < k < j, and 

3 

fc=l 

We also let r m = — x — v(a) and p(a) = ? — z^(a), where 



i/(a) 



if m is odd, 

6™ + i(a) if m is even. 



Now we are ready to introduce results of Sibuya [12] and Shin [TT] on existence of a solution 
that has a fixed asymptotic behavior in S-i U So U Si, and on its properties. 



Theorem 3. Equation @, wii/i a £ C m , Zias a solution f(z,a,X) with the following properties 

(i) f(z,a,X) is an entire function of z, a and X. 

d_ 

t)z 



(ii) f(z,a,X) and f'(z,a,X) = -§:f(z,a,X) admit the following asymptotic expansions. For 



each fixed e > 0, 

f(z, a, A) = z r ™ (1 + 0{z- 1 ' 2 )) exp [-F(z, a, A)] , 
f'(z, a, A) = -z r ™ + f (1 + 0{z- 1 ' 2 )) exp [-F(z, a, A)] , 

as z tends to infinity in the sector |argz| < ~~ £ > uniformly on each compact set of 
(a, X)-values. Here 



F(z,a,X) = +1 + Y - b 3 {a)z^ m+2 ^ 

v ' ' ' m + 2 ^ m + 2-2j < 3K ' 

i<i<f +i 



(Hi) For each 5 > 0, / and f also admit the asymptotic expansions, 

/(0, a, X) = [1 + O (A~')] A" 1 / 4 exp [L(a, A)] , (5) 
/'(0, a, A) = - [1 + O (A"")] A 1 / 4 exp [L(a, A)] , (6) 

as X oo in the sector | arg(A)| < tt — 5, uniformly on each compact set of a £ C m , where 



I mil j 

L(a,X) = £°° htm + P(t) + X-tf- ± bj (a)tf-J-^l 



dt. 



(iv) The entire functions X i— > /(0, a, A) and A \-$ f'(0, a, A) have orders p = A + — . 

Proof. See Theorems 6.1, 7.2, 19.1, and 20.1 in [12] for proof. Sibuya showed ([5|) and © with 
error terms o(l) and I improved them to O (X~ p ) in |11] , ■ 
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2.1 Spectral determinant 

In this subsection, we will introduce the function whose zeros are the eigenvalues of H. 
Let 



uj = exp 



2m 
m + 2 



and define 

G k (a) := (uj~ k a 1 ,oj- 2k a 2 , ■ ■ ■ ,uj- mk a m ) for k G Z. 

Then one can see from the definition of bj jk {a) that 

b jtk {G l {a)) = uj- j \ k (a) and bj {G e (a)) = u}~%(a) for £ G Z. (7) 

Also, by Theorem [3] there exists f (z , G k (a) , uj~ mk X) that decays in So and we let 

f k (z,a,X) :=f{u- k z,G k (a),u;- mk X). 

So /fc decays in S k and blows up in its two adjacent Stokes sectors. Since /o decays in Sq while f\ 
blows up in the same sector, we conclude that /o and f\ are linearly independent. Thus, any 
solution can be expressed as a linear combination of these two. In particular, 

f-i(z, a, A) = C(a, X)fo(z, a, A) + C(a, X)fi(z, a, A) for some C(a, A) and C(a, A). 

One can write these coefficients in terms of the Wronskians of solutions. That is, 

W-i,i(a,A) ~ W-i, (a,A) 

Vvo,i(o, A) Wo,i(a,A) 

where Wj t k = fjf'k~fjfk is the Wronskian of /j and f k . Since /o and /i are linearly independent, 
Wo,i(ct, A) ^ and likewise W_i i o(a, A) ^ 0. Then by Theorem E^i), for each fixed a G C m , the 
function A i— > C(a, A) is entire. Moreover, A is an eigenvalue of H if and only if C{a, A) = 0. We 
call C(a, A) the spectral determinant of H. 

2.2 Asymptotic expansions of C(a, A) 

In [10], the asymptotics of C(a, A) is provided and it is showed that all but finitely many zeros 
A of C(a, A) lie in a small sector containing the positive real A-axis. Here we will improve this 
asymptotics of C(a, A) as A —¥ oo in a sector containing the positive real A-axis as follows. 

Theorem 4. Suppose that m > 4. T/ien /or eac/i fixed a G C m and < 5 < ^^5, 

C(a,A)= [us + O (A" p ) ] exp [L(G _1 (a), w" 2 A) - L(a, A)] 

+ [wl+^W + O (A"') ] exp [L(G(a),a; 2 A) - L(a, A)] , (8) 
as A — > 00 in i/ie sector 

41 m I 4 

vr - + 5 < arg(A) < 7T ^— - 5. (9) 

If m = 3 iaen 

C(a,A) = [ - w~ 2 + O (\~ p ) } exp [L(G 4 (a), w~ 2 A) - L(a, A)] 

- [iws + O (\- p ) ] exp [-L(G 2 {a), co~ l A) - L(a, A)] , (10) 
as A —¥ 00 in i/ie sector 

+ 5 < arg(X) < tt - 5. (11) 

5 
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Proof. In Theorems 13 and 14 of Shin [10] . similar asymptotics of C(a, A) are proved with the 
error terms o(l) instead of O (X~ p ). With the improved asymptotics §5§ and (JSJ), one can closely 
follow proofs of Theorems 13 and 14 in jlO] to complete the proof. ■ 



Next, in order to further examine the improved asymptotics of C{a, A) in Theorem [H we will 
use the following asymptotics of L(a,X). 

Lemma 1. Let m > 3 and a 6 C m be fixed. Then 
L(a,X) =Y J K m j{a)\? +l ^ - ln(A) 

3=0 m 

b(- 1+—) 

as A — > oo in the sector | arg(A)| < tt — 5 with K m fi(a) = 2 2 ' / ™\ > and for j > 1, 

3 

K m j(a) = ^2 K m! j tk b j>k (a), (12) 
fc=Li^-J+i 

U m -I 1 



where for 0- < k < j 

J 111, — — u 



( 2 

m 



m,j,k 



if 3 = 1, 



( ln2 -Y^— 

\ ^ Is - 1 

\ 8=1 



2 /, „ nt^ 1 \ . m 

y 3 = tt + l> 771 ewen > 



m ^ 2s - 1 / 2 

\ s=l / 

1 „ /, ? — 1 ? — 1 1\ „ , , m 

—B[k-- , if j ^ l or j ^ — + 1. 

m V m m 2 / 2 



Proof. See the Appendix in |llj for a proof. ■ 

Remark 1. The coefficients K m j(a) have the following properties that can be deduced 
from (|12p and the definition of bj^{a): 

(i) The K m j(a) are all real polynomials in terms of the coefficients a of P. 

(ii) For each 1 < j < m, the polynomial K m j(a) depends only on a±, a%, ■ ■ ■ , aj. Furthermore, 
it is a non-constant linear function of a,- . 

3 Proof of Theorem [1] 

In this section, we will prove Theorem [TJ 

Proof of Theorem [TJ In order to get asymptotics for the eigenvalues A, we examine asymp- 
totics of C(a, A) since the zeros of C(a, A) are the eigenvalues of H. Since all but finitely many 
eigenvalues lie near the positive real A-axis, the asymptotics in Theorem [J] would be sufficient 
for our purpose. 

For m > 4 and a € C m fixed, we set C(a, A) = in ([8]), rearrange the asymptotic equation 
to get 

[1 + O {\~ p ) ] exp [L(G(a),uj 2 X) -L(G~ 1 (a),u~ 2 X)] = -u~ 2 ^ a \ 
and absorb the term [1 + (A~ p )] into the exponential term to obtain 
exp [L(G(a),uj 2 X) - L{G~ 1 (o),uj~ 2 X) + O (A^)] 
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or equivalently 

L(G{a),Lu 2 \) - L(G~ 1 (a),Lo^ 2 X) +0 (X~ p ) = (2n + l)m - ^^i, for some n € Z. (13) 



m + 2 



Since L(a, A) = -ff mj o( a )A p (l + o(l)) by Lemma[IJ one gets 



2vr\ 



L(G(o), w 2 A) - L(G _1 (o), w~ 2 A) + O (A _p ) = 2iK mfi (a) sin ^— J A'(l + o(l)), (14) 

as A —7- oo in the sector Q containing the positive real A-axis. Thus, the function 

A i — y H(a, A) := L{G(a),u 2 X) - L(G _1 (a), uT 2 A) + O (X~ p ) 

maps the region |A| > M and | arg(A)| < e for some large M > and small e > onto a region 
that contains the positive imaginary axis near infinity. Also, H(a, •) is analytic in the region 
|A| > M and |arg(A)| < e. Then, Sibuya [12j pp. 131-133] showed that for every closed ball 
{a € C m : \a\ < r}, there exists JVo = No(m,r) £ N such that for every integer n > Nq there 
exists exactly one A n for which (|13p holds, that is, 



(15) 



H(a,X n ) = 2n + l vri ^-1 

m + 2 

(Sibuya showed the existence of such an Nq with the estimate ([HI).) 
Hence, from Lemma Q] we infer that 

H(a, A re ) = L(G(a), u 2 X n ) - L{G~ l {a),uj~ 2 X n ) + O (X~ p ) 

m+l r ii- 
.2, 



K mJ (G(a))(a; 2 A n )^ - - K m , i (G- 1 (a))(a;- 2 A n ) 2+- 



E 

i=o 

1/ (G _1 (a)) „ " (G(o)) , / n. \ „ / . „\ 

+ — — In (ur 2 A.„) - v v " In (w 2 A n ) + O (X~ p ) 



I'll 



in 



m+l 



2i ^ Km,j(a) sin 



2(1- j)vr\ i+i-i , i/(a) 8ttz 



in 



X 



+ 



in m + 2 



(16) 



where we used v (G ±1 (a)) = ~v{a) and by © and for 1 < j < m + 1 , 



iT mj (G(a))(^ +1 ™ - ^(G-^a)) (^~ 2 ) 



2^ m 



tf mJ (a)c^'(a; 2 )^ +1 » - iiC mJ (a)^(^" 2 )^ + 1 " 



2iK m j (a) sin 



2(l-j> 



in 



(17) 



Next, one can combine (|15p and ()16p and rearrange the resulting equation to obtain ([3]). That is, 



m+l 



5^Cj-(o)A^ m +0(A n p )=n + - as n ^ 00, 

3=0 

where for < j < m + 1, 



Cj(a) := < 



-^mj(ffl) sin 



2v{a) 



m 



m 2 



(18) 



??? 



if j = — + 1, m even. 
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We still need to examine how many eigenvalues of H are not in the set {A n } n >Ar to complete 
our proof. To do this we will use Hurwitz's theorem in complex analysis. 

Hurwitz's theorem [1] says that if a sequence of analytic functions ip^ converges uniformly 
to an analytic function ip on any compact subsets of their common domain, then for all large k 
functions and (p have the same number of zeros in any open subset whose boundary does not 
contain any zeros of (p. Since the eigenvalues are the zeros of C(a, A) that is an analytic function 
in variables a and A, Hurwitz's theorem implies that eigenvalues varies continuously as a varies. 
That is, A n (a) varies continuously from A n (a) to A n (a) as a varies from a to 5 in the closed ball, 
where we used A n (a) in order to clearly indicate its dependence on a. 

When a = 0, there exists an integer N = N(m) < Nq such that the eigenvalues of H that are 
not in {A n (0)} n >Ar is exactly (iVo — N). Also, Hurwitz's theorem implies that when eigenvalues 
collide as a varies, the number of eigenvalues before and after the collision remains the same. So 
there is no sudden appearance or disappearance of the eigenvalues as a varies and hence, there 
are exactly (No — N) -eigenvalues that are not in the set {A n (a)} n >7v . Therefore, the eigenvalues 
{A n (a)} n >Ar satisfy ([3]). This completes the proof for m > 4. 

For m = 3, proof is very similar to the case m > 4 above and we will use fjlOf) in the place 
of ([8]). If C(a, A) = then we rearrange the asymptotic formula fjXOj) with C(a, A) = to obtain 

[l + O (X' p ) ] exp [-L(G 4 (a), u~ 2 \) - L(G 2 (a), w _1 A)] = e™, 

as A —> oo in the sector (1111). Next, like in (1171) we examine 



-K 3 ,(G\a)) (^ 2 A)^ + — -K 3 ,(G 2 (a)) (w _1 A)^ ~ 
= -[ w - 4 ^- 2 (5+V) + w " 2 ^-(^ +i 3 1 )]^3,i(a)A^ + ^ 

r • 2 (j-!) -i 1 i 1-3 

= -[e- J T i «-e- J T J »]if 3!3 ( fl )A2+T 



2isin v J ' K 3:j (a)\ 



where we used ([7]), fjl2[> . and ojz = —1. Then one can complete proof for m = 3 just like in the 
case of m > 4. ■ 

Remark 2. We order the eigenvalues in the order of their magnitude and we showed that 
eigenvalues vary continuously as the coefficients a = (a±, . . . , a m ) vary. For the large eigenvalues, 
A n (a) are continuous functions of a € C m (c.f. (|20p ). However, for small eigenvalues, even though 
eigenvalues varies continuously as a varies, we do not claim that A n (a) ordered in their order 
of magnitude are continuous. While these small eigenvalues are continuously moving in the 
complex plane, their indices can be switched and how we number these small eigenvalues does 
not affect the asymptotic result in Theorem [TJ 



4 Corollaries and proof of Theorem [2] 

In this section, we will introduce some results deduced from Theorem [TJ 

First, we invert (|3|), expressing A n as an asymptotic formula in n. Also, we obtain an asymp- 
totic formula for nearest neighbor spacing of eigenvalues and monotonicity of their magnitude. 



Corollary 1. One can compute numbers dj(a) explicitly such that 
m+i , s. % 

3=0 



m+1 / s. i» (l-J-) 

r , / 1 \ m + 2 V ml 4 ■ 

A n = 2^ dj(a) ■ I n+ - I +0(n ™+2) as n -> +oo. (19) 
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Also, the space between successive eigenvalues is 



2m n 
r-~ m — 2 

2/ \ m + 2 / \ rrr 

ml it \ [ 1 \ m + 2 



A n +i — A n = — T-j =-r- ■ I n H — ) + o(n m + 2 ) as n — > +oo. 

"+ 1 « n ^ +oom + 2^B(i ) l + i)y V 2; 1 > 

In particular, lim |A n +i — A n | = 00 and lim arg(A n ) = 0. Hence: 

n— >oo n—s>+oo 

|A n | < |A n +i| for all large n. (20) 
Proof. Since ([3]) is an asymptotic equation, one can solve it for A n to get (|19p . where do (a) = 

2m 

(cq) m + 2 and dj(a) for 1 < j < m + 1 can be computed recursively and explicitly. Then since 
do > 0, one can deduce 

lim arg(A n ) = and lim \X n \ =00. (21) 
Also, we obtain 



m — 2 



2771 2m / 1 \ m + 2 m — 2 

A n+ i = A n + —(do(a))^- (n + -) + o(n^) as n -> +00, (22) 

n->+oo m + 2 \ 1 J x ' 

where we use the generalized binomial expansion to (l + (n + ) a for some a G R. Since 
do > 0, equations (f2T|) and ([22]) imply limn^o |A n +i — A n | = 00 and hence, ([20]) holds. ■ 

When H is 'PT-symmetric, the eigenvalues either are real or else appear in complex con- 
jugates. Thus, if H is "PT-symmetric, then by monotonicity of magnitude of the large eigen- 
values ([20]) . H cannot have an infinite pairs of complex conjugate eigenvalues and hence all but 
finitely many eigenvalues are real. 

Next, we will prove Theorem Don a necessary and sufficient condition for an anharmonic 
oscillator H having infinitely many real eigenvalues. 

Proof of Theorem [2l Suppose that H or its translation is PT-symmetric. It is easy to check 
that the potentials V(z) and V(z — zo) generate the same set of eigenvalues for any zq G C. Thus, 
since H or its translation is PT-symmetric, eigenvalues of H either are real or come in complex 
conjugates. Then by (|20p . we conclude that all but finitely many eigenvalues are real. Since 
there are infinitely many eigenvalues [12] as we noted in the Introduction, there are infinitely 
many real eigenvalues. 

Suppose that H has infinitely many real eigenvalues. Then first, from the properties 
of K m j(a) in Remark [1] and (I18p . one can derive the following properties of Cj(a): 

(i) The Cj(a) are all real- valued polynomials in terms of the coefficients a of P. 

(ii) For 2 < j < m, the polynomial Cj(a) depends only on a\, a<2, . . . , a,j. Furthermore, it is 
a non-constant linear function of aj. 

Notice that if Cj(a) ^ K for some < j < m + 1, then ([3]) implies that all but finitely many 
eigenvalues are non-real. Thus, if H has infinitely many real eigenvalues, then Cj(a) & M for all 
< j < m + 1. 

Recall that co(a) is real and depends only on m, and ci(a) = 0. Next, if a\ 7^ 0, then we 
can choose zq = —^i and we consider the translation of H, replacing the potential V(z) by 
V(z — zo). The translated H has its potential with a vanishing z m_1 -term (i.e., a\ = 0). Then 
by the above properties (i) and (ii) of 02(a), 02(a) G M is a real polynomial in a\ and 02, and it 
is a nonconstant linear function in a<i- Thus, since a\ = G M and 02(a) G M, we obtain a2 G M. 
Likewise, since ai, 02, 03(a) G M, by the above properties of 03(a), 03 G M and hence, recursively 
we conclude that a G M m . Therefore, the translated H is PT-symmetric. ■ 
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